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Abstract 

We analyze solutions to a superconducting gap equation based on the two-dimensional 
Hubbard model with nearest and next-to-nearest neighbor hopping. The Cooper pair 
potential can be calculated exactly and expressed in terms of standard elliptic functions. 
The Fermi surfaces at finite temperature and chemical potential are also calculated based 
on the exact two-body S-matrix of the Hubbard model using the formalism we recently 
develop ed[7|], which allows variation of hole doping. The resulting solutions to the gap 
equation are strongly anisotropic, namely largest in the anti-nodal direction, and zero in 
the nodal directions of the Brillouin zone, but not precisely d-wave. For U/t = 13 and 
t'/t = —0.3, appropriate to BSCO, and a physically natural choice for the cut-off, our 
self-contained analytic calculations yield A an ti— nodaiA ~ 0-06 and maximum T c /t ~ 0.04 at 
optimal hole doping h = 0.15. For phenomenological fits to the Fermi surfaces for cuprates, 
we obtain the comparable value T c /t = 0.03 at optimal doping, both in good agreement 
with experiments. The superconducting gap is non-zero for all hole-doping h < 0.35 and 
increases all the way down to zero doping, suggesting that it evolves smoothly into the 
pseudogap. 
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I. INTRODUCTION 



The microscopic physics underlying high T c superconductivity in the cuprates is 
believed to be purely electronic in origin, and strongly correlated electron models 
such as the two-dimensional Hubbard model have been proposed to describe it l|. 
A partial list of more recent articles addressing the existence of superconductivity 
in the Hubbard model is 2|-|4j, and references therein. The Hubbard model simply 
describes electrons hopping on a square lattice subject to strong, local, coulombic 
repulsion. Since it is known that the condensed charge carriers have charge 2e, 
and thus some kind of Cooper pairing is involved, a central question has become 
"What provides the glue that pairs the electrons?" . This is especially puzzling since 
the underlying bare interactions are repulsive. The situation is completely different 
in the Bardeen-Cooper-Schrieffer (BCS) theory of ordinary superconductors, where 
the attractive glue is provided by the interaction of the electrons with the lattice 
phononsQ]. 

Since the Mott-insulating anti-ferromagnetic phase at half-filling is well under- 
stood, a large portion of the theoretical literature starts here and attempts to un- 
derstand how doping "melts" the anti-ferromagnetic order, and how the resulting 
state can become superconducting. This has proven to be quite challenging, perhaps 
in part due to the fact that anti-ferromagnetic order is spatial, whereas supercon- 
ducting order is in momentum space. Consequently this has led to many interesting 
and in cases exotic ideas, however the central question, "What is the glue?", and 
how it arises from strongly coupled physics, remains unclear. (For a review and 
other refereces, see [6|.) This suggests that it may be more fruitful to begin on the 
overdoped side, far away from any competing anti-ferromagnetic order, in order to 
understand the attractive mechanism in a pure form. Here the density is perhaps 
low enough that one can treat the model as a gas, with superconductivity arising as 
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a condensation of Cooper pairs as in the BCS theory, and we will adopt this point 
of view in the present work. The observation made in |t| now comes to bear on 
the problem: multi-loop quantum corrections to the scattering of Cooper pairs can 
actually lead to effectively attractive interactions, even though the bare model was 
defined with repulsive interactions. In the work the focus was on the thermo- 
dynamics at finite temperature and chemical potential, and evidence was presented 
for instabilities toward the formation of new phases as the temperature was lowered. 
However our purely thermodynamic formalism was unable to probe the nature of 
the ground states of these potentially new phases. The present work attempts to 
complete the picture. Namely, we explore how the attractive interactions described 
in {7] can lead to superconductivity, and what its basic properties are. 

Our starting point will be the BCS theory, but specialized to the Hubbard model. 
The original Cooper argument js] is quite robust, and shows that any attractive 
interactions near the Fermi surface lead to a pairing instability. We thus assume 
that the BCS construction of the ground state goes through, leading to the well- 
known gap equation^: 

A(k) = ~Iwf VM m^) tanh (£(k')/2T), £(k) = V((k) 2 + A(k)* 

(1) 

Here, A is the energy gap, E(h) is the energy of excitations above the ground state, 
and T is the temperature. We will later make some favorable checks on the approx- 
imations that lead to the above equation. In the above gap equation, V represents 
the residual interaction of Cooper pairs, and we will refer to it as the (Cooper) pair 
potential. The main new input is that we compute V from the Hubbard model, in- 
cluding quantum corrections, and show that it has attractive regions in the Brillouin 
zone. The other ingredient is £(k), which represents normal state quasi-particle ener- 
gies near the Fermi surface, where £ = at the Fermi surface. This can be identified 
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with the "pseudo-energy" in the thermodynamic approach described in These 
two ingredients lead to a self-contained analysis of the solutions of the above gap 
equation based entirely on analytic calculations carried out in the Hubbard model. 

An outline of the sequel, along with a summary of our results, goes as follows. In 
the next section we describe our conventions for the Hubbard model, with hopping 
strengths t,t' and the repulsive coupling U > 0, all with units of energy. The pair 
potential V is calculated in section III by summing multi-loop Feynman diagrams; 
the final result is expressed in terms of elliptic functions. It is demonstrated that, 
rather surprisingly, when U/t is large enough, there opens up a region of attractive 
interactions, i.e. negative V, near the half-filled Fermi surface. In section IV, the 



method developed in 
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9[ for the thermodynamics is reviewed, and Fermi surfaces 
are calculated as a function of doping for U/t = 13 and t'/t = —0.3, appropriate to 
the cuprate BSCO. The results in these sections III, IV constitute the main inputs 
for the study of the solutions of the gap equation, which is carried out in section V. 
The values we compute for the gap A and T c are in reasonably good agreement with 
experiments. The gap is anisotropic, in that it vanishes in the nodal directions and is 
largest in the anti-nodal, however it is not precisely of d-wave form. Our solutions to 
the gap equation persist to arbitrarily low doping, and we propose that they evolve 
into the so-called pseudogap, in accordance with recent experiments. In section VI, 
we repeat the analysis of the gap equation using a phenomenological fit to the Fermi 
surfaces. 

II. THE HUBBARD MODEL GAS 

The Hubbard model describes fermionic particles with spin, hopping between the 
sites of a square lattice, subject to strong local coulombic repulsion. The lattice 
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hamiltonian is 



H =~ t Yl «, a C rj , a )-t' Yl (fti,*^,*) +Uj2 n rt n ri (2) 

<ij>,a=t4 <i,j>' ,a=t4 r 

where r^r are sites of the lattice, < i,j > denotes nearest neighbors, n = c^c 
are densities, and c\c satisfy canonical anti-commutation relations. We have also 
included a next to nearest neighbor hopping term t', since it is not difficult to incor- 
porate into the formalism, and it is known to be non-zero for high T c materials. 
The free part of the hamiltonian, i.e. the hopping term, is easily diagonalized: 

#frcc = / d 2 k U k ^ C l,a C Ka ( 3 ) 
J a 

with the free 1-particle energy 

Wk = —It (cos(k x a) + cos(k y a)) — At' cos(fc x a) cos(k y a) (4) 

where t taken to be positive. In the sequel it is implicit that k is restricted to the 
first Brillouin zone, —it /a < k XjV < it /a, where a is the lattice spacing. 

Since the quartic interaction is local, we introduce the two continuum fields ip^^ 
and the action 

(5) 



S= f d 2 rdt(j2 i^l 9 ^* ~ H ) 



where H = H.f ree + 'Hint is the hamiltonian density. The field has the following 
expansion characteristic of a non-relativistic theory since it only involves annihilation 
operators, 

Mr) = j^ c ^ T (6) 

and satisfies at equal times {^ a (r), i^ a ,(r')} = 5(r—r')5 aj0l >. Since we have represented 
sums over lattice sites r as J d 2 r/a 2 , where a is the lattice spacing, c r = aip(r). The 
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interaction part of the hamiltonian is approximated as a continuum integral with 
density 

rtint = 2 nwjw ( 7 ) 

where tt = 2Ua 2 . Formally, the free part of the hamiltonian density is "Hf ree = 
S Q =t j,^!^ 7 r)^) where W is the differential operator W = u(k — > — iV r ), and is 
thus non-local. However this non-locality does not obstruct the solution of the model 
since the free term can be diagonalized exactly. The model can now be treated as 
a quantum fermionic gas, where the only effect of the lattice is in the free particle 
energies ui^. 

The field ip has dimensions of inverse length, and the coupling u units of energy • 

length 2 . In the sequel we will scale out the dependence on t and the lattice spacing 

a, and physical quantities will then depend on the dimensionless coupling 

u 2U 

Positive g corresponds to repulsive interactions. Henceforth, all energy scales, in 
particular, the single particle energies u^, the gap A, temperature, and chemical 
potential, will be implicitly in units of t. 

For both cuprates LSCO and BSCO, t » 0.3ev w 3000iT, U/t w 13, and t'/t 
approximately equals —0.1 and —0.3 respectively. Therefore, for most of the detailed 
analysis below, we set g = 26 and t'/t = —0.3 appropriate to BSCO. 



III. THE COOPER PAIR POTENTIAL V. 

The kernel V^(k, k') in the gap equation (CQ) represents the residual interaction of 
Cooper pairs of momenta (k, — k) and (k', — k'). It is related to the following matrix 
element of the interaction hamiltonian: 

V(k,k') = J d 2 v (k't,-k'||?{ int (r)|kt,-k|) (9) 
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To lowest order, V is momentum independent: V = g/2. 

In quantum field theory, the above matrix element of operators, in this case H^t, 
is generally referred to as a form-factor. Since there is no integration over time, 
this form-factor does not conserve energy, i.e. there is no overall 5-function equating 

to Ujfi. The form-factor can be calculated using Feynman diagrams as follows. 
More generally consider the form-factor (k 3 f, k 4 4- |"%int(r)|ki f, k 2 I). Represent 
the interaction vertex with two incoming arrows for the annihilation operator fields 
ip^^ and two outgoing arrows for the creation fields tf)l Furthermore, let such a 
vertex with a "node" • represent the operator "Hint- Then V is given by the sum over 
diagrams shown in Figure [U where p = (u, k) represents energy-momentum. 




P12 ~P P34-P P34-P 



FIG. 1: Multi-loop diagrams contributing to the Cooper pair potential V. 

Momentum is conserved at each vertex, however energy is not conserved at the 
vertex with a node. There is actually no fermionic minus sign associated with each 
loop since the arrows do not form a closed loop. Diagrams with a closed loop, such 
as in Figure El are zero because the integration over energy u inside the loop has 
poles in the integrand that are either both in the upper or lower half-plane, so that 
the contour can be closed at infinity without picking up residues. In other words, 
there is no "crossing-symmetry" as in relativistic theories. (For the contrary, see 
the non-zero loop integral C below.) This fact, which is unique to non- relativistic 
theories, allows us to calculate the kernel V exactly. At order g N+1 , specializing to 
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Cooper pairs ki = — k 2 = k and k3 = — = k', the diagram in Figure [T] factorizes 
and contributes 

(-^/2)^^ ^pll £(k)-£(kO Ar - (10) 
where £ is a 1-loop integral 

f dud 2 p ( % \ ( % \ . 

~ 7 (2tt) 3 \cu-co p + ie) \E 12 -cu -co kl2 ^ p + ie) 

where and k 12 are the total incoming energy and momentum, i.e. Ei 2 = + 
U]<_ 2 = 2wk and ki2 = ki + k 2 = 0. As usual, e is infinitestimally small and positive. 
The extra 1/2 N is due to the over-counting by allowing each loop to be £(k) or £(k'). 
Finally, summing over m, N gives 

^•^itWfw (12) 

It is important to note that the above V is exact and has a smooth g — > oo limit, i.e. 
it allows an expansion in the inverse coupling t/U, so is in a sense non-perturbative. 
One may be concerned that we formally summed a geometric series that potentially 
does not converge. In answer to this, there are certainly regions where £ is small 
enough that the series converges. Also, this summation is known to give the correct, 
exact, S-matrix for non-relativisitic quantum gases, and this S-matrix has all of the 



right properties in the strongly coupled unitary limit [10], namely, it gives the correct 
diverging scattering length at the renormalization group fixed point, and the bound 
state. The only difference here is that the the kinetic energy k 2 /2m is replaced with 
cuk for the Hubbard model, which does not affect these arguments. 

The w-integral can be performed by deforming the contour to infinity, giving 
£(k) — i J <i 2 k/87r 2 (u;k — u p +ie). Note C is imaginary as e — > 0, thus in the formula 
(fl2"|) . C is really the imaginary part of C as e — > such that V is real. Then, integral 



over p can be performed analytically 



£(k) = - 



(w k - 4f ) 



ix VK + 4 + 4t')K(4 - Wk) + i6f (f - 1; 



1/2 K ( 16(^-1) 



(w k + 4t') 2 - 16 



(13) 



with Wk Wk + ie, where K is the complete elliptic integral of the first kind. Note 
that the momentum dependence of the kernel only enters through the variables Wk, 
i.e. V'(k,k') = V(w k ,w k ,). 




FIG. 2: Diagrams with closed loops, i.e. arrows circulating in the same direction, vanish. 

Non-zero solutions to the gap equation possibly signifying superconductivity can 
only arise if the effective interactions are attractive, i.e. if the kernel V is negative. 
For g small and positive, the effective coupling V remains repulsive. However, as 
pointed out in [7|, for g large enough, V can become negative in certain regions 
of the Brillouin zone. Since we are interested in a small band of energies near the 
Fermi surface, V(u,u') with u' = u is a suitable probe of these attractive regions. 
In Figure [3] we plot this V for g = 10 and 26 at fixed if = —0.3. One sees that for the 
smaller g, V is everywhere positive, however for larger g it flips sign. As explained in 
more detail in [3] , this feature is reminiscent of what occurs near the fixed point of 
quantum gases in the unitary limit, where for the same analytic reasons, the effective 
interactions can be either attractive or repulsive depending on which side of the fixed 



point of the BEC/BCS crossover lOj. Using the formula f fT3|) . one can show that there 



is a region of negative V for g > 13.2. This minimal value of g depends on if and this 
dependence was studied in [7| based on the formula ( TT3"j) . Around g = 13 — 15, the 
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attractive region is a narrow band[7(|. It will also be instructive to view a contour 
plot of V in the first Brillouin zone, see Figure HJ 
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FIG. 3: Plot of the Cooper pair potential V(uj,uj) for g = 10,26 and t' = —0.3. (Colored 
photos on-line). 

The main effect of a non-zero t' is the following. For g large enough, V(uj,u) 
becomes negative for u > At'. Thus when t' is negative and |t'| large, attractive 
interactions exist deeper inside the half-filled Fermi surface. If superconductivity 
indeed arises from these attractive interactions, then a non-zero if can play a signifi- 
cant role, otherwise the attractive interactions only exist too close to the vicinity of 
the half-filled Fermi surface where it has to compete with the known Mott-insulator 
phase. There is actually some evidence that superconductivity does not exist for 
f = 
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FIG. 4: Contour plot of the Cooper pair potential V(uJk, Wk) for 9 = 26 and t' = —0.3 in the 
first Brillouin zone, with axes k x ,k y . (All subsequent contour plots in the Brillouin zone 
follow the same conventions.) In the central light region, the interactions are repulsive, 
whereas in the colored regions attractive. Values of V can be inferred from Figure[3j (Color 
online) . 

IV. THE FERMI SURFACES AS A FUNCTION OF DOPING 



We will utilize the approach to the thermodynamics of particles developed in 




which is based on a self-consistent re-summation of the exact 2-body scattering. The 
occupation numbers are parameterized by two pseudo-energies e-y,(k), which satisfy 
2 coupled integral equations with a kernel related to the scattering of spin up with 
spin down. For equal chemical potentials, due to the SU(2) symmetry, = e±, and 
both occupation numbers are given by 

/« = ^ihr < 14 > 
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where e satisfies the single integral equation: 

f d 2 k' 1 
E(k)=u k -fi-J ^C(k,k / ) e£(k;)/T + 1 (15) 

The kernel G is related to the logarithm of the 2-body S-matrix, and is built 
from the same ingredients as the kernel V in the gap equation, since it also involves 
a sum of Feynman diagrams of the kind shown in Figure [TJ It is somewhat more 
complicated than the pair potential V since the total incoming momentum is not 
zero. Namely, consider the same loop integral as in the previous section but with 
k x + k 2 ^ 0: 

L(k ll k 2 )=i/ 7 ^ — = X + z 7 (16) 

where X and 7 are defined to be real. Then the kernel takes the following form 

g < bgf'/M-g/f) (17) 



21 °\l/g R + a/2, 
where the renormalized coupling is = g/{l — g r y/2). (We are not displaying 
the momentum dependence; it is implicit that G = G(k = k^k' = k 2 ).) The 
renormalized coupling is related to the gap equation kernel of the last section as 
follows: gn(k, — k) = 2V"(k, k). The quantity X represents the phase space available 
for 2-body scattering. The argument of the log is the exact 2-body S-matrix. 

We define the hole doping h as the number of holes per plaquette, which is related 
to the density n as follows: 

7 (27p e^/ T + l ( 8) 

where a is the lattice spacing. 

The integral equation (fT5|) was solved numerically using an iterative procedure, as 
explained in [7 1 . The solution for the pseudo-energy yields the relation between the 
chemical potential \i and the hole doping h. For a given h, fi(h) of course depends 
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h 




FIG. 5: Hole doping h as a function of chemical potential fi at the reference temperature 
T = 0.1. (g = 26, t' = -0.3.) 

on temperature, but only weakly j^]. For the subsequent analysis we determine fi(h) 
at the low reference temperature T = 0.1. The result is shown in Figure [51 

As described in jjj], for low enough T, there are regions of fi,T where there are 
no solutions to the integral equation at low enough T. Regions of the non-existence 
of solutions are very similar to that shown in Figure 8 in[7 1 , For hole dopings < 
h < 0.25, there are no solutions for temperatures T c in the range 0.02 < T c < 0.07. 
(This is why we chose the reference temperature To = 0.1 to be above these potential 
transition temperatures.) It was suggested in [3] that the non-existence of solutions 
could represent an instability toward the formation of a new phase, however the 
nature of these new phases cannot be determined based on what we have done 
so far; one needs a complementary bottom up approach that is based on the zero 
temperature ground state. This is the subject of the next section. As we will see, 
the critical temperatures inferred from this thermodynamic analysis are consistent 
with the critical tempertures computed from the gap equation in the next section. 
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The Fermi surfaces are the contours <Sf(/-0 that are solutions to e(k, /i) = 0, 
where \x depends on hole doping as in Figure [5j The pseudo-energy e thus represents 
the quasi-particle dispersion relation in the normal state. These Fermi surfaces are 
shown in Figure [6] for < h < 0.4. These calculated Fermi surfaces are in reasonably 
good agreement with experiments. Comparison with Figure [TU1 which is based on 
a phenomenological fit to the data, suggests that one needs to include additional 
hopping terms in the bare hamiltonian, such as next-to-next neighbor. We will refer 
to wave- vectors k that point to (tt, tt) (and 90° rotations thereof) as being in the 
nodal direction, whereas those pointing to (0, n) as in the anti-nodal direction. 

In the same Figure |6] we also display the region of attractive interactions based on 
the Cooper pair potential V calculated in the last section. Before even solving the 
gap equation, one can make some predictions concerning the existence of supercon- 
ductivity based on the attractive region of V. Namely, for high enough hole doping, 
about h > 0.3, there is no attractive region near the Fermi surface, and thus no 
superconductivity. It is also clear from Figure [6] that the regions of the Fermi surface 
in the anti-nodal direction are the most important since this is the direction with the 
greatest overlap with the attractive region. As we will show in the next section, this 
feature is primarily responsible for the anisotropy of the gap A(k), and explains why 
the gap is zero in the nodal direction, at least for moderately high doping h > 0.1. 

V. SOLUTIONS TO THE GAP EQUATION 

The gap A(k) is defined and measured for k along the Fermi surface. However 
near half-filling, there is no k from the origin that intersects the Fermi surface in the 
anti- nodal direction; see Figure El Consequently, polar plots of A(k) for k originating 
from the center of the Brillouin zone, and covering all of it, are potentially misleading, 
since A(k) is not well-defined in the anti- nodal directions. In fact, this also implies 
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FIG. 6: Calculated Fermi surfaces for hole doping < h < 0.4 in steps of 0.1, with h = 0.4 
the innermost curve. The grey (pink on-line) region corresponds to the attractive region 
for g = 26, t' = —0.3 based on the Cooper pair potential V. 

that A cannot be strictly d-wave as denned from the origin, for example, it cannot 
be of the simple form A oc | cosk x — cosk y \, where k is measured from the origin. 
Thus, in the first quadrant of the Brillouin zone, it is more convenient to work with 
the vector k originating from the node (it, it), i.e. k = k„- )7r + k, where k^ = (it, it). 

In the integral gap equation (JT]), one must integrate over a narrow band around 
the Fermi surface Sp(fi). Let 9 be the angle of k relative to the horizontal line 
through the (it, it) node: k = —K,(cos(8)i + sin(#) j). If the integration is over a 
narrow band 5Sf of width 2 5k around the Fermi surface, then the gap equation 
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takes the following form in the first quadrant: 

m Jo Jk f {6')-5k &{K,V) \ 11 J 

(19) 

where E = +\f£, 2 + A 2 , and Kp(6) is the length of /c along the Fermi surface. Here 
V_ is the negative (attractive) part of V since only when V is negative are there 
solutions; the repulsive parts of V are already incorporated in determining the Fermi 
surfaces. The doping h dependence of the above equation is implicit in kf(9). 

It remains to determine the cut-off 5k. There is some arbitariness in the choice 
of 5k, and this is the weakest aspect of our calculation since T c certainly depends on 
it, just as T c depends on the Debye frequency in ordinary superconductors. It could 
be viewed as a free parameter that needs to be fit to the data. Or, one could carry 
out a sophisticated renormalization group analysis requiring A to be independent of 
5k, but this is at the expense of introducing an arbitrary scale that has to be fit to 
experiments. Instead, we found the following choice to be physically meaningful and 
well-motivated. As in the BCS theory, 5k should be related to the properties of the 
potential V itself, namely the region in which it is attractive relative to the Fermi 



surface 



181 ] . Define Ky such that V is attractive for k < Ky in the nodal direction. 



For g = 26, t' = —0.3, Ky w 2.7. We then take 5k as a measure of the distance of 
the Fermi surface to the edge of the region of attractive interactions in the nodal 
direction: 5 k = \Ky — kf(tt /4)\. From the Figure El one sees that 5 k is quite small, so 
that our choice does correspond to a narrow band. For the Fermi surfaces computed 
in the last section, 5k w 0.04 for doping h — 0.15 and this value will be used in the 
subsequent analysis. 

The function £(k) in the gap equation represents the quasi-particle dispersion 
relation of the normal state, and is zero along the Fermi surface. We wish to carry 
out a self-contained calculation, thus we equate £ with the pseudo-energy e of the 
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last section. On the other hand, one can infer £ from experiments, and we will repeat 
the analysis with such a £ in the next section. 

A measure of the validity of the BCS approximation is the combination of the 
effective coupling V and the phase space integration in the gap equation, i.e. the 
parameter v = —k,p5kV/87t. From Figure [31 V ~ —20, and also Kp m 2.7, 5k ~ 0.04. 
This gives v ~ 0.05, which appears to be sufficiently small for the BCS approximation 
to be valid. 

We solved the gap equation numerically by discretizing the integrals and solving 
the resulting set of coupled non-linear equations with Mathematica. In Figure [7] 
we plot the zero temperature gap as a function of the Fermi surface angle 9 for 3 
different hole dopings. One sees that A(k) is highly anisotropic in the Brillouin 
zone: it is largest in the anti-nodal directions, and zero in the nodal direction for 
high enough doping. At low doping, the anisotropy is less pronounced. This property 
can be traced to the detailed shape of the Fermi surface in comparison to the region 
of attractive interactions displayed in Figure El More specifically, for high enough 
doping, the Fermi surface does not intersect the region of negative V in the nodal 
direction, however it always does in the anti-nodal direction. This effect is even more 
pronounced for the experimentally determined Fermi surfaces, as will be described 
in the next section. 

Our solutions to the gap equation are not exactly d-wave, more specifically, are not 
proportional to | cos k x — cos K y \, which is approximately Kf \ cos 29\/2, for the Fermi 
surface approximated as a circle of radius Kp. Experimental data indicates a gap 
closer to the d-wave form, however some data does show a tendency for it to flatten 
out in the nodal direction, as in our solutions. However the precise shape of the gap as 
a function of 9 will change if the cut-off 5k is made to depend on 9, instead of a simple 



constant as we have done here. Some features in Figure [7] are reflected in the data 14] 



in particular, the central region around 9 = 45° where the gap is smallest widens 
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with increasing doping. We also wish to point out that our Figure |6] is suggestive of 
an observation made in [l^]: the Bogoliubov quasiparticle interference, indicative of 
the existence of Cooper pairs, disappears alon g th e diagonal line connecting the two 
anti-nodes (0, 7r) and (vr,0). (See Figure 3 in 14j.) Interestingly, this diagonal line 
is very close to the contour that separates attractive from repulsive regions of the 
pair potential, however we are unable to make a direct connection at present. These 
observations are more pronounced in Figures [11] and [TUJ 
A 
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FIG. 7: The zero temperature gap A as a function of the Fermi surface angle 6 for hole 
doping h = 0.10,0.15,0.20. 

The gap in the anti-nodal direction is plotted as a function of temperature for 
h — 0.15 in Figure El Where it goes to zero defines T c , in this case T c « 0.04. Our 
results for the zero temperature gap in the anti-nodal direction and T c for various 
doping are summarized in the table below. Figure |9] plots both the zero temperature 
gap and T c as a function of doping. Moving down from the overdoped side, the 
maximum T c occurs first at h — 0.15, in good agreement with experiments. For 
t = 3000-fT, one obtains A = 19meV and T c = 120/i at h = 0.15, compared with 
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FIG. 8: The gap A in the anti-nodal direction as a function of T for hole doping h = 0.15. 
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As explained in the last section, there is a succinct reason for why there is no 
superconductivity at high enough doping, roughly h > 0.35, since beyond this, no 
part of the Fermi surface overlaps with V_ . See Figure [61 A slightly lower value of 
h > 0.3 is more typical in experiments. However there is no mechanism in the gap 
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FIG. 9: The zero temperature gap A in the anti-nodal direction and T c , in units of t, as a 
function of doping h. 

equation to turn off the gap at low enough doping, and it continues to increase all 
the way down to zero doping. Superconductivity would turn off if g were smaller, 
namely around g = 15, since in this case the region of attractive interactions is a 
narrow band, and at low enough doping the Fermi surface does not overlap with 
it [3| - However if eq. (JS]) is accurate, g is roughly twice as large, so this does not 
account for the disappearance of superconductivity at low doping. Although this may 
appear problematic when one compares with the usual phase diagram of the cuprates, 
there is growing experimental evidence that this is actually the correct behavior 1^, 



151 ] . Namely, it has recently been found that the d-wave superconductivity gap 



evolves smoothly into a d-wave pseudogap whose magnitude continues to increase to 
arbitrarily low doping, where superconductivity is not present. In other words, the 
so-called pseudo-gap energy scale T* may be the continuation of the superconducting 
gap, had there been no other other mechanisms to destroy it. As stated explicitly in 



15] . these results are inconsistent with 2-gap scenarios. One check of this is that if 
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one identifies T* with A, then T* = 195K for t = 3000-fT at doping h = 0.05, again in 
reasonable agreement with experiments. This suggests that on the underdoped side, 
superconductivity is perhaps destroyed by competition with other orders, presumably 
ant i- ferromagnetic, bringing T c to zero, even though the gap A is still physi cally 



present. It could also be destroyed by phase decoherence, as suggested in[16|, [17]. 
These effects are of course not implemented in our gap equation, and it is beyond 
the scope of this paper to address this, for example by comparing the free energies 
for the competing orders or to study to phase fluctuations. 

VI. SOLUTIONS TO THE GAP EQUATION FOR PHENOMENOLOGI- 
CALLY DETERMINED FERMI SURFACES 

In this section we repeat the analysis of solutions of the gap equation with the 
normal state quasi-particle dispersion relation £(k) determined from experimental 
data, but with the same pair potential V computed in section III. There is extensive 
data on the Fermi surfaces for the compound E^S^CaC^Os+s. A tight-binding fit 
to the data was performed in [3] based on the data in [? ]. The result is that the 
Fermi surfaces are the contours £(k) = for the following function: 

£(k) = — /2 — 2 (cos k x + cos k y ) + 0.6513 cos k x cos k y — 0.4455 (cos 2k x + cos 2k y ) 
— 0.1716 (cos 2k x cos k y + cos k x cos 2k y ) + 0.6357 cos 2k x cos 2k y (20) 



We have rescaled the result for £ in 12j so that t — 1. The parameter ju serves as 
a renormalized chemical potential. For hole dopings h = 0.07,0.08, 0.14,0.17,0.19, 
$ = -0.0500, -0.178, -0.503, -0.688, -0.809 respectively. 

Let us assume that the underlying Hubbard hamiltonian has the same U/t = 13, 
i.e. g = 26, and still only has nearest and next-nearest neighbor hopping parameters 
t and t', where from eq. ( 1201) one reads off t' = —0.163. The other terms in eq. (I20jl 
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should be viewed as generated by the interactions, for example by equations such as 
in section IV for the pseudo-energy e. In Figure [10] we display the resulting Fermi 
surfaces against the region of attractive interactions V- as computed in section III, 
but with t' = —0.163. In comparison with Figure El one sees that the effect that 
leads to the anisotropy of the gap is more pronounced: the Fermi surfaces are pulled 
away from the attractive region in the nodal direction in a stronger manner, which 
implies that the gap will continue to be zero in the nodal direction for lower values 
of h in comparison to the last section. 




-3-2-10 1 2 3 



FIG. 10: Tight-binding fit to the Fermi surfaces for Bi2Sr2CaCu20g+5 at dopings 
0.07,0.08,0.14,0.17,0.19, eq. (f20|) . The grey (pink online) region corresponds to the at- 
tractive region for g = 26, t' = —0.163 based on the Cooper pair potential V. 
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In Figure [TT] we plot the solution to the gap equation for h = 0.07, 0.14. The same 
cut-off 5k = 0.04 as in the last section was used. For the reasons stated above, the 
gap is zero for a wider region centered at 8 = 45°. Plots of the g cip clS db function of 
temperature are very similar to those of the last section, and lead to a slightly lower 
T c , i.e. T c 0.03 at optimal doping h = 0.14. 
A 




9 (degrees) 



FIG. 11: The zero temperature gap A as a function of 9 for hole doping h = 0.07,0.14 
based on the phenomenological £ eq. ([20]) . 



VII. CONCLUSIONS 

The central proposal of this work is that quantum loop corrections to the Cooper 
pair potential, as computed here in the two-dimensional Hubbard model, could be 
responsible for the effectively attractive interactions near the Fermi surface that 
lead to the phenomenon of high T c superconductivity. The validity of this idea is 
easily explored, since the pair potential can be calculated exactly, and the results 
favor our proposal. We showed that the resulting analysis of the solutions to the 
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superconducting gap equation leads to definite predictions for the anisotropy of the 
gap, its magnitude, and T c , all in reasonably good agreement with experiments. 
We explained in a clear manner why superconductivity disappears at high enough 
doping; in our calculation h > 0.35. 

We found that the non-zero solutions to the gap equation continue undiminished 
on the underdoped side, all the way down to zero doping, and this appears to be 
consistent with recent experimental results QQ, with the interpretation that the 
d-wave superconductivity gap evolves smoothly to the d-wave pseudogap. If the 
attractive interactions considered in this paper are indeed responsible for supercon- 
ductivity, then this feature suggests two items on the more speculative side: 

• On the underdoped side, superconductivity is perhaps destroyed by competition 
with the anti-ferromagnetic phase which is known to exist at very low doping. Here 
one must bear in mind that we kept only the attractive part of the pair potential for 
the superconducting gap equation, but inside the Fermi surface the interactions are 
still largely repulsive. Another possibility is that it is destroyed by phase decoherence 
of the gap HQ- 

• The pseudo-gap energy scale T* may thus represent the hypothetical contin- 
uation of the superconducting gap had it not been destroyed by the mechanisms 
suggested above. If one identifies T* with A, then T* = 210K for t = 3000fT at zero 
doping, again in reasonable agreement with experiments. 

If these ideas are correct, then the emphasis should shift from trying to under- 
stand "doping the Mott insulator" to its opposite, that is to say, understanding 
how populating the superconducting state can destroy it due to the competing anti- 
ferromagnetic order, phase decoherence, or perhaps something else. This issue has 



been studied experimentally in significant detail[14|. The latter approach may be 
more tractable if based on the concrete description of high T c superconductivity 
presented in this paper. 
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